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Abstract 

We study the 3d effective field theory of a weakly coupled two Higgs-doublet scalar 
model at high temperature. Our model has three scalar quartic couplings and an 
0(4) symmetry which is spontaneously broken by a nonzero scalar field vev at zero 
temperature. Using high temperature perturbation theory, renormalization group 
arguments in 4 — e dimensions, as well as 1/N expansion techniques in 3 dimensions, 
we argue that the transition from the high temperature symmetry restoring phase to 
the low temperature phase is first order for a range of scalar couplings. This result 
is not due to gauge couplings. We discuss the implications of our results for the 
transition in the two Higgs-doublet electroweak model, especially when the lightest 
neutral Higgs is heavier than the W-bosons. 
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Introduction. The topic of high temperature phase transitions in 4d scalar models 
has attracted a lot of interest recently due to the possibility of baryogenesis during the 
weak phase transition. It has been argued that baryogenesis requires a sufficiently strong 
first order phase transition [1,2,3,4]. 

The nature of the phase transition can be determined if one has a reliable determination 
of the finite temperature effective scalar potential V e ff. f\ It has been known for some time 
[7,8,9] that at sufficiently high temperature the ground state of a 4 model is symmetry 
preserving with < <fi >= 0, even if the zero temperature ground state is not. What is 
not well understood is how the transition from the high temperature phase to the low 
temperature phase proceeds in generic models with several couplings. 

One can imagine two possibilities as the temperature is lowered from a very high value: 
i) the minimum remains at < >= until a temperature T\ at which point it jumps 
discontinuously to nonzero values | < <p > I 7^ 0, or ii) the minimum remains at < 
<f) >= until a temperature T 2 at which it moves continuously away from the origin as 
the temperature is lowered further. Case i) indicates a first order transition and case ii) 
indicates a second order transition. 

The determination of the transition order is not always easy due to infrared divergences 
in high T field theory which can render perturbation theory useless in calculating V e ff too 
close to the origin and at temperatures near the transition temperature. In this case, one 
can i) try resuming perturbation theory or ii) use renormalization group (RG) methods 
or iii) use 1/N expansion techniques to gain some insight. When a model has several 
couplings, any one method cannot be trusted to give a reliable description of the transition 
for all values of the (perturbative) couplings. 

The most thoroughly investigated electroweak model, the minimal standard model, can 
be argued [10-14] to give a first order transition for Mm ggs sufficiently smaller than M\y. 
In this case, high T perturbation theory can give a good description of the transition if 
used carefully. 

The leading expression for the high temperature effective potential when Mmggs <S Mw 

^What we mean by V e ff is the sum of one-particle irreducible diagrams with no external legs. This 
quantity need be neither real nor convex [5] , in contrast to the standard definition of the effective potential 
given in many textbooks [6]. However, the minimum of the real part of our effective potential still gives 
the ground state of the theory. 
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is of the form [7,15,10-14] 

V; //= m 2 (T)|$| 2 -5T|$ 3 | + ^|$| 4 , (1) 

where $ is the Higgs doublet, 5 is proportional to the cube of the SU(2) gauge coupling 
and m 2 {T) oc (T 2 — T 2 2 ). This potential describes a first order transition whose minimum 
jumps discontinuously to a value oc 5T/X at a temperature 7\ > T 2 . RG methods combined 
with e-expansion techniques [15,16] as well as computations of V e ff using 1/N techniques 
(when My/ ~ NMniggs) [17,18] are consistent with this result. 

With the current experimental bounds on Mm ggs [19], the situation Mm ggs > O(Myy) is 
the one of physical interest. Unfortunately, for Mmggs > Mw, perturbation theory breaks 
down and there is really no clear picture of what happens. However, one expects that at 
best the phase transition will be very weakly first order. 

It has been suggested [11] that one should augment the minimal standard model with 
extra scalars in order to produce a big enough first order effect for values of the Higgs mass 
larger than the current experimental bound. The idea is to enlarge the ratio S/X during 
the phase transition. 

One suggestion [11], is to add an electroweak singlet scalar S with a large enough 
coupling to the Higgs doublet <3> of the form £$t$S' 2 . It has been argued that if S is 
massive enough it contributes to V e ff(T = 0) in a way that reduces the value of A without 
decreasing M Higgs . 

Another nonminimal model is the two Higgs-doublet electroweak model. Due to the 
number of unknown parameters in this model it makes sense to study a toy model before 
attempting to understand the nature of the transition in the general case. In this letter we 
will study such a toy model. In particular, we study the 4d weakly coupled scalar model 
with tree potential 
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Vfa, fa) = - v 2 ) 2 + (0 2 - v 2 ) 2 + ^-0 2 2 + ^ • 2 ) 2 , (2) 



£ ,2 J.2 , a 
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where A > £ (to ensure real physical masses) and <f>i and 2 are both 0(4) vectors. Unless 
otherwise stated, we assume all couplings are arbitrarily small (their ratios, however, may 
not be) and nonnegative. 

The model (2) has an 0(4) symmetry as well as the discrete symmetry 0i <-> 02- It 
is also invariant under 4>i ~^ ~4>i or 02 - ► —02- The potential (2) is a subset of the 



general potential describing the two Higgs-doublet model which has five quartic couplings 
[20], not three. One of the additional couplings breaks the Z 2 symmetry (f) 1 <-> <f> 2 , the 
other is allowed because the potential needs to be invariant only under SU(2) x U(l), not 
0(4) ~ SU{2) x SU{2). 

Although the 2 Higgs-doublet model has been studied at one-loop by previous authors 
[21], there has never been a reliable demonstration of a first order transition when the 
lightest neutral Higgs weighs at least 0{Mw)- 

In this article we argue that perturbation theory gives a reliable description of a first 
order transition if a/ A and cn/£ are sufficiently big. We present additional arguments for a 
first order transition based on RG group methods combined with the e-expansion, as well 
as computation of V e ff based on 1/N expansion techniques. We stress that the first order 
behavior is not due to gauge fields. We will later address the implications of our results 
for baryogenesis in the two Higgs-doublet electroweak model. 

We study not the full 4d model, but the 3d effective field theory describing the high T 
limit of the 4d model. Before proceeding with our analysis we introduce all the tools we 
will use in the context of a simple 4 model. 

Effective 3d Theory. Let us consider the ungauged scalar theory with 4d potential 

VX0) = ^(0 2 -. 2 ) 2 - (3) 

The model is, at high T, formally equivalent to a euclidean field theory with one compact 
dimension. For physics at scales less than 0(2nT) it is sufficient to study the effective 3d 
Lagrangian that results from integrating out the compact dimension. The 4d fields can be 
expanded as 

<j>(x,T)=J2Mx)r(r), (4) 

n 

where r parameterizes the compact dimension and the ^"(O) = ^"-(T -1 ). Thus the 4d 
field yields an infinite tower of 3d fields when the compact dimension is integrated out in 
the action. Furthermore, for what we are interested in, only the zero model n — 0, for 
which iP°(t) = ip°(0) is important in the effective 3d model because the others pick up 
large 0(2irT) masses from the 4d kinetic term. 

The largest effect of integrating the nonzero modes out from the theory is to give an 
0(AT 2 ) correction to the mass-squared of the zero mode [16]. Therefore, we study an 
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effective 3d theory with potential 

Vb) = ±{T-A&)<f + ~<p\ (5) 
where we have defined the 3d quantities 

^ = 00/^, A = AT, v 2 = v 2 /T. (6) 

The effective theory can be studied by computing the 3d V e ff, or by RG methods. 
The ungauged model is believed to have a second order transition [22]. This result is not 
accessible by means of high temperature perturbation theory. To see this we need only 
compute V e ff to one-loop. 

The field dependent mass for the potential (5) is 

m 2 (^T) = A(i2^ + r-4iJ 2 ). (7) 

We expect the 3d model to give a good description of V e ff(ip,T) for m 2 < Att 2 T 2 , i.e. for 
T 2 near Av 2 and § V < 4vr 2 T 2 . Q 

The one-loop contribution to the effective potential is then (m 2 > 0) 

l m , r ^9 9t Am 2 Iml 3 , , 

-Tr ln[-<9 2 + m 2 ] = - + const. (8) 

2 ' Air Ylix 

We have used a sharp momentum cutoff to regulate the integral. 

The linearly divergent field dependent term is absorbed by renormalizing the 3d theory. 
The field dependent part of the one-loop effective potential is then 

Veff = —(T - 45 V - — + V- (9) 
e}} 48 V ^ 12tt AV v ' 

This result reproduces the leading temperature dependent part of the so-called daisy sum in 

the high temperature field theory. Taken at face value, (9) indicates a first order transition 

with a discontinuous jump in the minimum | < > | oc \f\T at a temperature T\ slightly 

above T 2 = 2v. Since we have already indicated that the phase transition is believed to 

be second order, eq. (9) must not give a reliable determination of V e ff. In fact, it is well 

*We will freely switch back and forth from 3d quantities and 4d quantities, dropping the zero mode 
subscripts. 
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known [8,9] that higher loop corrections are significant. In 3d, the scalar coupling A has 
dimension 1. The dimensionless perturbation parameter is actually ~ A/m, a result that 
can be checked by examining multiloop graphs [8,9]. At T = T2, A/m ~ VXT/cJ), and 
for the apparent minimum at | < <f) > | oc v^T this expansion parameter is of 0(1) and 
therefore perturbation theory cannot be relied upon to distinguish between 1st and 2nd 
order behaviour. 

In the case of an 0(N) symmetric model described by (3), V e ff has been computed 
using 1/N expansion techniques [23,24,22,8] to next-to-leading order. The results indicate 
that the phase transition is second order [24,17]. This method quickly becomes technically 
complicated, so that V e ff has not been computed beyond next-to- leading order (known 
results should not be trusted to give a reliable description of what happens for low values 
of N). 

Fortunately, there have been many studies of scalar models based on the renormalization 
group method [22]. A trajectory in coupling constant space which leads to an infrared stable 
fixed point under the RG is identified with 2nd order behaviour. There are two possibilities 
to study the RG flow equations for a 3d theory. One is to work in 4-e dimensions to some 
order in e. Although one is interested in e = 1, even lowest order results in e can give 
useful insight into the nature of the phase transition.^ The other is to stay in 3d and to 
work to some order in 1/N. As with 1/N computations of V e ff, leading order results in 
this expansion are generally only expected to give a good description of the model when 
there are many scalars. 

For a simple 0(N) symmetric scalar theory, RG equations have been studied to a high 
order in both the 1/N expansion and the e-expansion [22]. Both results indicate a 2nd 
order transition and, in both cases, leading order results give a faithful indication of this. 
For example, for N = 1 and 4d potential (3), the 0(e) and 0(A 2 ) RG flow equation is 
[16,22,26] 



eA -T7TT A2 - ( 10 ) 



dX 

~dt 16tt 2 

Here, t increases in the infrared. This has an infrared unstable fixed point at A* = and 
an infrared stable fixed point at A* = 167r 2 e/3. In the infrared limit, every positive A is 



§ The absence of an infrared stable fixed point at 0(e) although not a proof [25] that the model has 
first order behaviour, may hint at such behaviour for some range of couplings. 
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driven to this point. 

For an abelian Higgs model, i.e. a gauged 0(2) model, with gauge coupling g and 4d 
potential (3) where 2 = • 0, the RG equations to 0(e) have no infrared stable fixed 
point [15,16], suggesting that a first order transition may be possible. In fact, high T 
perturbation theory can be used to give a reliable description of a 1st order transition for 
g 2 ^> A ^> g 4 [15,7]. Perturbation theory is reliable because the perturbative expansion 
parameters go like some (positive) power of X/g 2 at both minima at T — T\ [10,12,14]. 

For Mm g g S > Mw the above analysis breaks down. For the abelian Higgs model, a com- 
putation of V e ff to next-to-leading order [17,18] in the 1/N expansion gives a second order 
transition. Futhermore, calculations of critical exponents during a 2nd order transition to 
next-to-leading order in 1/N remain physical for values of iV > 10 [29,15] suggesting that 
the prediction of a second order transition for M w < M Higgs and N > 10 is reliable. In 
fact, the phase transition is known to be second order for iV = 2, M w < M Higgs due to 
other reasons [30]. This demonstrates that the RG analysis to 0(e) in 4 — e dimensions 
can fail to indicate 2nd order behaviour. 

As an aside, we note that much less is well understood when the gauge group is non- 
abelian. As with the abelian case, there is no infrared stable fixed point to 0(e) in 4 — e d. 
Assuming that infrared divergences in loops due to the self interactions of the magnetic 
fields are cutoff by a nonperturbatively generated T-dependent mass [27] , it is again pos- 
sible to argue [10,12,13] that for g 2 3> A the transition is first order. For g 2 < O(A), 1/N 
studies [18,31] suggest that the transition may be second order for a model with enough 
scalars. For a model with just a few scalars, recent results [32] due to resuming high T 
perturbation theory in the minimal electroweak model also suggest 2nd order behaviour 
for g 2 < O(A). In addition, it has been argued [33] that the existence of a fixed point at 
0(e) in 2 + e d may give a reliable indication of a second order transition at d — 3 and low 
N. 

The Toy Model. We are finally ready to demonstrate that the ungauged model 
described by (2) has a 1st order phase transition for a range of couplings. To our knowledge, 
what we present below is new. 

To begin, we examine the RG equations for (2) in 4 — e dimensions. To O(e), one has 
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[26] 

dX , 1 



rfa 1 ,, 3 2 , / , 

— — = ea -\Aa + Ea-\ — a k (11, 



This has the following fixed points: 

i) a* = 0,C = 0,A* = 0, 



ii) a* = 0,C = 0,A* = 4vr^e, 

94 48 

m) a* = 0, f = — vr 2 e, A* = — n 2 e. (12) 
) ' ^ 13 13 K J 

Linear stability analysis reveals the following facts. The fixed point i), the Gaussian 
fixed point, is unstable in all directions. The fixed points ii) and Hi) are both unstable in 
the a direction: in the infrared limit, small (positive) values of a are driven larger. The 
lack of a stable fixed point suggests that the purely scalar model may undergo a first order 
transition. Indeed, the fact that the coupling a is responsible for the lack of stable fixed 
points suggests that any first order behaviour will be due to a. 

We therefore examine high T perturbation theory for a ^> A and a ^> £. To simplify 
the analysis, we consider 

a = / 3 , e = / 3 , « = / 2 , (13) 

for small positive /.Q This choice of couplings also ensures [14] that all mass scales of 
interest are much less than T, thus allowing us to work with an effective 3d theory. Under 
such circumstances, we may now argue that perturbation theory gives a reliable description 
of a first order transition. The argument follows closely the one for the abelian Higgs model. 
The effective 3d potential that we study is 

V tree (<Pn <P2) = + Vlf + ^f(<Pi ■ ^) 2 + + <pl)M*, (14) 



• Our particular assignments are only for numerical convenience and we will later make clearer for what 
values of couplings we predict first order behaviour. 
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where 

M 2 = ^^T°-^ = £{l + 0(f))(T«-T°), (15) 
Tl = Q\v 2 /(a + 5\) = Qfv 2 (l-0(f)), (16) 

and ipi, tp2 are the (properly normalized) zero modes of the 4d fields 0i, 02- 

As with the simple scalar model discussed previously, we expect the effective 3d theory 
to give a reliable description of the phase transition as long as all mass scales are less than 
2ttT. 

The field dependent mass matrix has the eigenvalues 

fp + + 0(f), f P - + 0(f), fcr + + 0(f), fa. + 0(f), (17) 



where 



p± = ^(T 2 -T 2 )±±<j )l .(f> 2 , 

a± = —(T 2 -T 2 ) + - 
± 36 V 2J 6 



^ + ^± v /(0? + 0i) 2 + i2(0 1 -0 2 ) 2 



(18) 



The first two eigenvalues each have multiplicity 3. 

The one-loop contribution includes a linearly divergent correction which renormalizes 
the 3d theory. The remaining finite contribution is 

V^oo P = (34 + 3p f + 4 + a!) + 0(f). (19) 

At one-loop, we have V e ff = Vt ree + V\-i oap . At T = T2, the real part of the potential has 
a minimum at | < <pi > \ = \ < 02 > | ~ yg and < 4>i ■ 4>2 >— 0(f). 

Taken at face value, to lowest nontrivial order in /, V e ff implies that the minimum is 
away from the origin by the time the mass at the origin vanishes and hence that a first order 
transition has occurred. The dimensionless loop expansion parameters are proportional to 

>*V„£1 and (20 ) 
m m m mm 

Here m is a field dependent mass eigenvalue. At T = T 2 , all mass eigenvalues are zero 
at the origin so that V e ff is not to be trusted near the origin. To verify that (23) gives a 
reliable prediction of a first order transition we must first demonstrate: i) there is another 
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minimum degenerate with the one at the origin at some T = Ti > T 2 , and ii) the loop 
expansion parameters are small at both minima at T — 7\. 

Since the 0i • 2 terms are 0(f 2 ) at the minimum at T = T 2 , they do not affect the 
lowest order (in /) minimization conditions for < <f>i > and < 02 >• It turns out that for 
the purposes of discussing the leading order (in /) characteristics of the phase transition 
one is allowed to just ignore these terms. To lowest order in /, the 3d effective potential 
we study for T close to T 2 is therefore 



(PT), 2 2 , 2 M \ 2 2n / 

Veff = ^TT^l + vl) + -IT&l + <P%) 



3 r ' ■ " 3 



4| vti 1-2/ 2 vt-1 12?r 



^(T 2 -T|) + i(0? + ^) 



(21) 



Although this contains a cubic term only at T = T 2 , it still leads to first order behaviour 
because for T close enough to T 2 the last term behaves like a cubic term oc (<p\ + 4>\)^- 
One can verify that Ti = T 2 + 0(f)T 2 and | < 0i > | = | < 02 > | = O(T) at the second 
minimum at T = Ti. At the origin, all mass eigenvalues go as f 3 T 2 2 and thus the loop 
expansion parameters go as fy/J and yff. Away from the origin, in the 0i-02 = direction, 
the mass eigenvalues only increase and the expansion parameters decrease. Therefore, for 
sufficiently small /, the perturbative expansion is well controlled at both minima at T = Ti. 
Just below T = Ti the global minimum is away from the origin; this is the signal of first 
order behaviour. 

We end our analysis of the model (2) by discussing three issues: 

1) There are values of 0i • 2 for which some of the field dependent mass eigenval- 
ues of (14) vanish and hence perturbation theory cannot be reliably used to compute 
V e ff((f>i, 2 , 0i • 02) for such values. For example at T = T2, six of the 0(f 2 ) mass eigen- 
values (17), (18) vanish at 0i • 02 = suggesting that perturbation theory is not reliable at 
the apparent minimum away from the origin. Actually, at T = T 2 and 0i • 02 = these 
mass eigenvalues go like 0(/ 3 T 2 2 ) so that perturbation theory is reliable at the minimum. 
However, there are still values 0i • 02 = 0{f(f>\ , /0|) for which the mass eigenvalues vanish 
and perturbation theory must be resummed in order to get a more reliable V e ff. 

2) So far we have said nothing about the transition for A ~ £ ~ a. In this case, since 
perturbation theory breaks down, we must resort to other methods. Even though the RG 
analysis in 4 — e d has no infrared stable fixed point to 0(e) the transition may be second 
order for these values of the couplings. To learn more we can use 1/N expansion techniques 
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[22,23] when 0i,02 are A^-vectors. We have studied two interesting limits: 



i) £ = A and (iV 2 A), (Net) fixed for increasing N, 
ii) a = £ = A and (NX) fixed for increasing N. 



(22) 



At large N, case (22z) leads to first order behaviour while case (22ii) leads to second order 
behaviour. 

To show this last statement it is convenient [23,17,18] to perform the rescalings 0i —> 
a/]V0i, 02 — ► VN4>2 and v — > V^/Vv. Then, the 3d effective theory is described by the tree 
potentials 



i) V 



S (A+a) -iV 



ii) V = N 

+0(1 



T\ A 2 
72 A -12 V 



(23) 



for A — > \/N 2 ,a — > a/A" and A — > A/A/", respectively. The 3d kinetic terms for both cases 
are also proportional to N, due to the rescalings of the fields. 
To compute V e ff in the 1/N expansion one performs the shift 



<pi -> vi + (pi/y/N, (p 2 -> ^2 + w/Vn, 



(24) 



dropping all terms linear in the quantum fields 0i,0 2 and computes consistently all one- 
particle-irreducible diagrams with zero external quantum fields to some order in 1/AT. The 
computations are simplified if dimension two auxiliary fields are first used to eliminate the 
quartic couplings [22,23]. 

Since the scaling (22i) only confirms our perturbative analysis we will just state how 
the analysis proceeds. V e ff at O(N) depends only on the combination 0i ■ 2 both at tree 
level and after all the O(N) diagrams are included. The real part of V e ff is an increasing 
function of 0i • 02 and therefore in the large A" limit the vacuum energy is minimized at 
0i ' 02 — 0. In this limit, 0(1) corrections are easy enough to compute; they include the 



cubic term — 

LZ7Y 



ft) 



which leads to first order behaviour. Note that unlike the 



perturbative result (21) the cubic term contains no T 2 
term is lower order in N. 



T$ term. This is because such a 
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To study the scaling (22ii) we add the following terms involving dimension two auxiliary 
fields, x+ an d X-i to the tree potential (23ii): 

3N ( X+ + X- AT 2 \ 2 3N ( X+ - X- AT \ 2 

-^(^1 + ^2)-^ -TTT^ o 5-^1- (25) 



2 AT V 2 g vri r^/ y 2AT V 2 3 

after which the potential becomes to O(N) 

3N , „ 9x JV , N 3JV, , /l, 9 9x M 2 \ , x 

= ~iAT (x + + X ] + T^ 1 ' " X ] + ~T (x+ + X ] [ 6 (V? ' + ^) + AT J ' ( 26 ) 

where M 2 = T 2 A/36 — \v 2 /6. After the shift (24), the Lagrangian is only quadratic in 
the quantum scalar fields; these can then be integrated out at one-loop to give the leading 
order effective potential in terms of y?i,y?2 and x±- Eliminating the auxiliary fields from 
this effective potential by their equations of motion gives V e ff(<pi, ^2) to O(N). 

The eigenvalues of the scalar field mass matrix of (26) are just x+ an d X— Both 
eigenvalues have multiplicity N. Therefore, the O(N) quantum corrections are 

^Tr \n[-& + xA + f Tr ln[-^ + X -\- (27) 

The field dependent divergent part can be absorbed by renormalization of the tree potential 
and the remaining field dependent finite part yields ( X ± > 0) 

Veff = V - — (xl + X 2 - ) • (28) 
Here, V is given by (26). The equations of motion for x± are then 

X ± = ^(0i ± 2 ) 2 + M 2 — ^-VXI. (29) 
6 LZtc 

It is now straightforward to verify several facts. For M 2 > 0, X ± > and the minimum 
of the potential is at < 0i >=< 2 >= 0. For M 2 < 0, x± are both zero at 4>\ ■ 2 = 
and A(0 2 + 2 )/6 + M 2 = 0. Now since = d{ ^ 2 y oc X ±, the point X+ = X- = 

has the interpretation [22,23] as the symmetry breaking minimum for M 2 < 0. This point 
moves continuously away from the origin as T 2 falls continuously below 6v 2 ; this is the 
signal of a second order transition. 

Although we have only demonstrated that (2) has a second order transition for the 
particular case A = £ = a and at large N, we might expect a similar outcome for N = 4 
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and a < 0(A, £), i.e. the first order behavior driven by a is wiped out by the quantum 
fluctuations due to A, £ for small enough a. A completely reliable determination may only 
be possible on the lattice. 

3) Finally, in this section, we determine more precisely the ratio a/ A for which pertur- 

2 

bation theory breaks down completely and the maximum ratio \ for which the effective 



3d theory can be trusted. We examine the one-loop 3d potential for <f>i ■ 
dropping for the time being the T 2 — T| part of the "cubic" term: 



V, 



eff 



(XT), 2, 2N2 , M % 2 , 2n a * 



<t>l) 



0, A = £ and 



(30) 



M 2 and T 2 2 are given by (15) and (16), respectively. For this potential, T 2 — T| = 
a 3 T 2 /9A7r 2 (a + 5A), and M(Xi) = ^a 3 /ATi/187r . 

Now note that the high T perturbative expansion parameters for (2) are only propor- 
tional to \T/m,^T/m,aT/m. To get an estimate of the constants of proportionality we 
examined higher loop corrections in the 3d theory which grow with the number of scalars. 
For £ = A, the large N loop expansion parameters for N = 4 are AT/(67rm), aT/(67rm), 
where m 2 is a field dependent mass eigenvalue. At 7\ these expansion parameters are ~ 1 
for a ~ 5A. In addition, at T — T\ and for a > 5A the T 2 — T| term is less than the <f)f + (p 2 . 
term in the last term of (30), i.e. (30) gives a good description of the potential. 

The condition that all field dependent mass eigenvalues are less than 2irT at the mini- 
mum at T = T 2 gives a 2 < 127r 2 A. Thus, we predict first order behaviour for around" 

..2 

(31) 



a , a 
T2^ <A< 5' 



For £ = this last lower bound on a is instead around A < a/3. 

The Two Higgs Doublet Model. As we have mentioned, the general two Higgs 
doublet model has 5 scalar quartic couplings, not three. In addition, it has two scalar vevs 
v i, v 2 , gauge couplings, and dimension two terms which softly break the discrete symmetry 
0i — ► —0i- The scalar potential is typically written [20] 



V = Ai($l$i -vlf + A 2 ($ 2 $2 



.2\2 



+ A 3 + $2$2 - V l ~ V 2 



+ X 4 



1 1 We remark that as long as all couplings are perturbative at zero T these conditions for first order 
behaviour are not unnatural. 
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+A 5 Re(${$ 2 ) - videos C +A 6 Im(<l>{$2) - v x v 2 sinC 



(32) 



where the Xi,( are real parameters and $i,$2 are complex Higgs doublets. We have the 
following identifications (when Ai = A 2 , A 4 = A 6 = 0) 



between (2) and (32) with the normalization = 2 /2, $2$ 2 = 0|/ 2 - 

It is reasonable to ask what can be learnt about the the transition in the two doublet 
model from the analysis of our toy model. One can convince oneself that (32) will display 
first order behaviour with the following assumptions: i) v\ = t>2, ii) A5 is big enough 
compared to the other quartic scalar couplings and gauge couplings, and iii) as with the one 
doublet model, there is a nonperturbatively generated mass for the nonabelian magnetic 
gauge fields [27]. Conditions i) and ii) ensure that the mechanism for first order behaviour 
(within perturbation theory) is the same as for our toy model and condition iii) is necessary 
to ensure that the loop expansion parameter associated with the nonabelian gauge fields 
is small at small scalar field values at T = T x . Although our toy model (2) has <fri ■ $2 = 
at tree level, the fact that the 0(e) RG results are independent of terms quadratic in the 
fields in the tree potential strongly suggests (but does not prove) that we need not require 
cos£ to be small for (32) to exhibit first order behaviour (cos^ 7^ 1 generically implies CP 
violation in the Higgs sector). The perturbative analysis for cos( ~ 1 is more subtle than 
for our toy model. 

However, note that for sin( pa 0, A6(Im$|$ 2 ) 2 — ^A6(0iT0 2 ) 2 , with nonzero entries 
T 12 = T 3 4 = — T 2j i = — T4 3 = 1. We have checked that this term can play the same 
role as the last term in (2). Thus, we predict first order behaviour when sin( pa and A6 
(which determines the pseudoscalar mass Ma) is sufficiently larger than the other couplings, 
even when the mass of the lightest neutral Higgs is near or above Myy. The precise 
determination of the parameter space for which first order behaviour is possible will require 
a more careful analysis, however the ratio M^/M/j could be as low as 3. 

Acknowledgements. We are grateful to Rob Leigh and especially Raman Sundrum 
for interesting and useful discussions. 

Note Added. We have just received an interesting article on the electroweak phase 
transition with a singlet: J.R. Espinosa and M. Quiros, IEM-FT-67/93, [hep-ph/9301285| . 



Ax + A 3 = A/6, A 3 = £/6, A 5 = 2a/3, 



(33) 



13 



References 

[1] M. Dine, P. Huet, R. Singleton and L. Susskind, Phys. Lett. B257: 351 (1991). 

[2] M. E. Shaposhnikov, Nucl. Phys. B287: 757 (1987); 

A. I. Bochkarev and M. E. Shaposhnikov, Mod. Phys. Lett. A6: 417 (1987); 

L. McLerran, M. Shaposhnikov, N. Turok and M. Voloshin, Phys. Lett. B256: 451 (1991); 

N. Turok and J. Zadrozny, Phys. Rev. Lett. 65: 2331 (1990). 

A. Cohen, D. Kaplan and A. Nelson, Phys. Lett. 245B: 561 (1990) and Nucl. Phys. B349: 
727 (1991). 

[3] M. E. Shaposhnikov, Phys. Lett. B277: 324 (1992) and Erratum, ibid B282: 483 (1992). 

[4] M. E. Shaposhnikov, JETP Lett. 44: 465 (1986), Nucl. Phys. B287: 757 (1987) and Nucl. 
Phys. B299: 797 (1988). 

[5] G. Slade, Comm. Math. Phys. 102: 425 (1985) and 104: 573 (1986); 
A. Dannenberg, Phys. Lett. B202: 110 (1988). 

[6] C. Itzykson and J.-B. Zuber, Quantum Field Theory, McGraw Hill (1980). 

[7] D. A. Kirzhnits and A. D. Linde, Phys. Lett. B42: 471 (1972) and Ann. Phys. 101: 195 
(1976). 

[8] L. Dolan and R. Jackiw, Phys. Rev. D9: 3320 (1974). 

[9] S. Weinberg, Phys. Rev. D9: 3357 (1974). 

[10] M. Dine, R. G. Leigh, P. Huet, A. Linde and D. Linde, Phys. Rev. D46: 550 (1992) and 
Phys. Lett. B283: 319 (1992). 

[11] G. W. Anderson and L. J. Hall, Phys. Rev. D45: 2685 (1992). 

[12] C. G. Boyd, D. E. Brahm and S. D.H. Hsu, CALT-68-1795, HUTP-92-A027, EFL-92-22. 

[13] M. E. Carrington, Phys. Rev. D45: 2933 (1992). 

[14] P. Arnold, Phys. Rev. D46: 2628 (1992). 

[15] B.I Halperin, T.C. Lubensky and S. K. Ma, Phys. Rev. Lett. 32: 292 (1974). 

[16] P. Ginsparg, Nucl. Phys. B170: 388 (1980). 

14 



[17] V. Jain, MPI-Ph/92-72 and Jiep-ph/9209227 



[18] V. Jain and A. Papadopoulos, LBL-33067 and |hep-ph/921027S| , to appear in Phys. Lett. B. 



[19] ALEPH, DELPHI, L3 and OPAL collaborations, as presented by T. Mori, ICHEP, August 
8, 1992, Dallas Texas. 

[20] J. F. Gunion, H. E. Haber, G. Kane and S. Dawson, The Higgs Hunter's Guide, Addison- 
Wesley (1990). 

[21] A. I. Bockharev, S.V. Kuzmin and M.E. Shaposhnikov, Phys. Lett. 244B: 275 (1990); 
N. Turok and J. Zadrozny, Nucl. Phys. B369: 729 (1992). 

[22] J. Zinn-Justin, Quantum Field Theory and Critical Phenomena, Oxford University Press 
(1989) and references therein. 

[23] S. Coleman, R. Jackiw and H. D. Politzer, Phys. Rev. D10: 2491 (1974). 

[24] R. G. Root, Phys. Rev. D10: 3322 (1974). 

[25] T. P. Cheng, E. Eichten and C. F. Li, Phys. Rev D9: 2259 (1974). 

[26] J.H. Chen, T.C. Lubensky and D.R. Nelson, Phys. Rev. B17: 4274 (1978), 

[27] A. D. Linde, Rep. Prog. Phys. 42: 389 (1979); 

T. Appelquist and R. D. Pirarski, Phys. Rev. D23: 2305 (1981). 

[28] W. Buchmiiller, T. Helbig and D. Walliser, DESY 92-151, Nov. 1992. 

[29] S-K. Ma, Phys. Rev. D7: 2172 (1973). 

[30] C. Dasgupta and B. I. Halperin, Phys. Rev. Lett. 47: 1556 (1981). 
[31] S. Hikami, Prog. Theor. Phys. 64: 1425 (1980). 

[32] J. R. Espinosa, M. Quiros and F. Zwirner, CERN-TH.6577/92, December (1992). 
[33] J. March-Russell, Phys. Lett. B296: 364 (1992). 



15 



